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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

A COMPARISON OF THREE THEORETICAL METHODS OF CALCULATING 
SPAN LOAD DISTRIBUTION ON SWEPT WHIGS 
By Nicholas H. Van Dorn and John DeYoung 


SUMMARY 

Three methods for calculating span load, distribution, those 
developed by V.M. Falkner, Wm. Mutterperl, and J. We is singe i’ have 
been applied to five swept wings . The angles of sweep ranged from 
—4-5° to +45°. These methods wore examined to establish their 
relative accuracy and ease of application. Experimentally determined 
loadings were used as a basis for judging accuracy. For the 
convenience of the readers the computing forms and all information 
requisite to their application are included in appendixes. 

* From the analysis it was found that the Woissingor method 

would bo bost suited to an over-all study of the effects of plan 
form on the span loading and associated characteristics of wings. 

The method gave good, but not best, accuracy and involved by far 
the least computing effort. The Falkner method gave the best 
accuracy but at a considerable exponso in computing effort and 
honce appeared to bo most useful for a detailed study of a spocific 
wing. The Mutterperl method off ©rod no advantages in accuracy or 
facility over oithor of the other methods and honce is not recommended 
for use . 


INTRODUCTION 

In an effort to roach higher fli^it speeds, designers are 
turning to widely divorsifiod types of plan forms the aerodynamic 
characteristics of which are as yet unkn own . Since the multiplicity 
of such designs precludos an experimental investigation of occh, 
considerable attention has boon directed toward moans of obtaining those 
characteristics theoretically. Usually the basis for such theoretical 
investigations is span, loading. While the precise span loading 
itself may not be considered of major importance, it is believed 


2 


NA.CA EM No. A7C31 


that any method giving reasonably accurate predictions of span 
loading would bo amenable to simple extensions which would givo 
reasonably accurato values of such, characteristics as lift-curve 
slope, spanwiso center of pressure position, downwash at arbitrary 
locations, and rolling moments due to sideslip or rolling. 

A number of methods have boon developed for predicting the 
span loading of swept wings of arbitrary tapor and aspect ratio, 
but very few attompts have been mado to compare, for sovoral methods, 
predicted and experimentally measured loadings on identical wings. 

The investigation reported heroin was undertaken to provide such a 
comparison of prodicted and measured span loadings. The theoretical 
methods have been evaluated in terms of rolativo accuracy, manner 
and consistency of orror, and tediousnose of application. 

The mothods developed by V.M. Falkner (reference l), Wm. Muttorporl 
(reference 2), and J. Woissinger (reference 3) have been applied to 
five wings produced by swooping the wing panels of an airplane through 
a rango of -45° to +4-5° (reference 4) . The span load distributions 
so calculated have been compared with those obtained experimentally 
at the time of the investigation of reference 4. In addition, tho 
lift— curve slopes and spanwise contor-of— prossure position prodictod 
for oach wing by the several methods have been ccmparod witn those 
values obtained experimentally. 

Throughout the calculations a chock was made of tho timo 
required for oach method and for tho various ports of each method. 

From those observations a comparison was -made of tho rolativo 
todiousness of oach method, and indications wore obtainod as to 
which parts might be rendered less difficult and time consuming. 

Finally, in order to enable immediate application of the methods, 
all necessary tables, computation forms, and step-by-step computa- 
tion instructions for each are included in the appendixes. It is 
believed that with these aids a computing staff could undertake the 
computation of swept- wing characteristics with little additional 
supervision. In addition, for the convenience of the reader, there 
are included In the appendixes any mathematical derivations or 
developments not immediately obtainable from the references . 
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SYMBOLS 

General 

wing area, square feet 
effective 1 wing span, feet 
effective aspect ratio (t> 2 /s) 
semispan (b/2), feet 
wing chord, feet 
root chord, feet 
average chord (s/b), feet 

taper ratio, tip chord divided by root chord (ct/c 0 ) 

sweep angle of quarter-chord line positive for sweepback, 
degrees 

geometric angle of attack of wing measured from angle for 
zero lift, degrees 

geometric angle of attack of wing root section, degrees 

local geometric angle of attack, degrees 

longitudinal coordinate of downwash point positive 
forward, feet 

lateral coordinate of downwash point positive to 
right, feet 

dimensionless lateral coordinate of downwash 
point (y/a) 


-•-Ih all instances except the unswept wing, the actual nip chord was 
not parallel to the wind stream. An effective tip chord that 
was parallel was therefore assumed such that Hie wing area 
remained constant. The effective span is the span to this 
effective tip. 




MCA EM No. A7C31 


longitudinal coordinate of vortex element positive 
forward, feot 

lateral coordinate of vortex element positive to 
right, feet 

dimensionless lateral coordinate of vortex element (y/s) 
chord at epanwlse station q, feot 
density of air, slugs per cubic foot 
air-stream velocity, foot per second 

air-stream dynamic pressure (§pT 2 ) , pounds per square foot 

lift, pounds 

lift coefficient (L/qS) 

section lift coefficient 

vorticity, foot per second 

circulation, feet squared per second 

spanwise center of pressure position 

differential pressure between upper and lower surfaces 
of wing, pounds per square foot 

static pressure, potrnds por square foot 

free— stream static pressure, pounds per square foot 

pressure coeff icient[(p~p 0 )/q ] 

induced vortical velocity or downwash, feet per second 

downwasli anglo, the ratio of downwash to free-stroam 
velocity 

spanwise position in circular coordinates (cos^q 
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Symbols Pertaining to the Falkner Method 

y v semispan of horsoshoo vortox (s/20), foot 

x*' dimensionless longitudinal coordinato (x/yv) 

x* dimonsionloss longitudinal coordinate of control point 

relative to vortex 

x’ longitudinal coordinate referred to 0.5c lino, foot 

8 circular longitudinal coordinate (cos”* 1 2x r /c) 

y*‘ dimensionless lateral coordinato (y/yy) 

j* dimonsionloss lateral coordinato of control point relative 

to vortox 


a mj n unknowns in distribution series 
M number of vortices in chordvise direction 


V designates which of M vortices in chordvise direction 


A,B,C functions used in development 


P v,A; 

r v,B; 

r v> 


c irculation increment of vortox in two-dimensional flow, 
feet squared por second 


Gv,A; 

G v,B; 

G v,C 


dimensionless circulation factors C^v.a/AVj^v.b/ISVj 

r v,c/CV) 


r v total circulation of vortex in two-dimensional flow, 

feet squared per second (?v,A + ^v,B + ^v,c) 

H dimensionless lateral coordinate of midpoint of 

specific vortox (y/s) 

circulation of specific vortex in three-dimensional 
flow, foet squared, per second 


* 
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Ar.tn swoop anglo of loading edge, positive, for swoopback, 
degrees 


Symbols Portaining to the Muttorporl Mothod 

b* dimensionloss span along tho 0.25c lino (b/c 0 cos A) 

s' dimensionless semispan along the 0.25c line (b*/2) 

y* dimensionless coordinate of control point along 

line parallel to 0.25c line (y/c 0 cos A) 

y T dimensionless coordinate of vortex element along 

0.250c line (y/c 0 cos A) 

Br, B L perpendicular distance from. 0 . 25 c line to control point 
divided by c Q 

Ar, Al distance along 0.25c line from center section to base 

of perpendicular to control point divided by c 0 «' 


5 Ar - y l 

^ cos -1 y’/s 1 



Symbols Pertaining to the Weissinger Method 
ar local aspect ratio (b/c) 

G(<p) dimensionless circulation (P(y)/bV) a continuous 

function of rf 
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m 


M 

n 

V 

V- 

P-1 

OtL 

Cpv 


number of stations at which speoific circulation 
is to "be determined and at which downwash is 
summed 

number of stations at -which f n ^ and lA(v,p.) 
is to he determined J ‘ 

denotes at which of m points specific circulation 
ordinate occurs 

denotes at which of m points dovnwad3 is summed 

denotes at which of M points f n y_ or or 

La(v,{-U ordinate occurs 

denotes which of m terms in interpolation function 
f n,[i 

circular coordinate of point n (rnt/m+1 * cos -1 t}) 

circular coordinate of point v (va/m+1 = cos -1 t)) 

circular coordinate of point p (pjr/M+1 = 003 “ 1 rj) 


Gn 

Gy 

c v 

ar v 


} 


dimensionless circulation at spanwise station Cfc 

dimensionless circulation at spanwise station 
<Pn = <Pv 

chord at spanwise station 9 y 
specific local. aspect ratio (b/c v ) 

influence function 


La(ti^) 

L a (v,p) 

®v,n jB*v,ni B# v > 

^Vjn^v v 

b*yb*v |i'gv f -^ unc ^ oriS us ®^ in mathematical development 

SavJSv 
t/T(t);K 


f n,p.^n,p 


Pv 


interpolation function used in mathematical solution 
l/ary 
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DESCRIPTION OF WINGS INVESTIGATED 

The five wings to which the methods have teen applied are 
described in reference 4. Briefly, they were produced by sweeping 
the wing panels from an existing airplane to five angles of sweep. 
Each wing then consisted of a center section, the two main panels, 
and the two tip sections . The airfoil sections of the root and 
tip were generated by direct extension of the surface of the panels. 
The geometric characteristics of the wings are as follows: 


r A 

X 

AR 

(effective) 

s 

-4-5.2° 

0.376 

2.99 

335.5 ft 2 

-29.6° 

.405 

4.45 

282.3 ft 2 

•9° 

.542 

4.47 

201.8 ft 2 

31.0° 

.442 

4.66 

288.4 ft 2 

46.4° 

.418 

3.45 

309.5 ft 2 


In all applications presented herein It was assumed that all 
section lift-curve slopes were 0.103 per degree, the average value 
of this parameter for the sections at the ends of the unswept wing 
panel. Actually, the local section slope varied from root to tip; 
however, because of the nature of the sections generated by 
extending the wing panels, exact values of this function could not 
be determined. Corrections to the theoretical methods to account 
for such a variation were omitted from the computations, although 
the effects of such an emission are discussed later. 

Ibe loading on a wing can be separated Into the basic loading 
(that oxisting at zero over-all lift) which is a function of twiBt, 
camber, flap deflection, and plan form; and additional loading, 
which is a function of plan form and angle of attack. For purposos 
of analysis in this report, attention has boon directed solely 
towards the additional loading. The wings experimentally investigated 


2 To agroo with tho definition of sweep used in the theoretical 
methods of span loading prediction, sweep haB been referred to 
the sweep of the line joining the quarter— chord points at root 
and tip. In reference 4 tho sweep was defined somewhat 
differently . 
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wore essentially devoid of any camber or twist. Any evidence of 
basic loading shown by experiment was removed from the loading 
curvos used as a basis of comparison. Thus, for purposes of analysis, 
the ving was rcplacod by a flat plate and local angles of attack 
becomo synonymous witli over-ell angles of attack. The character- 
istics of wings having cambor or twist tho variation of which is 
free from discontinuities, however, could be determined equally 
well by any of the methods simply by using the true local angle 
of attack (as measured from tho angle of zero lift) at each point 
considered. Further discussion of this problem is given in the 
appondixes . 


PROCEDURES 

All methods described heroin are extensions of simplified 
wing theory and so are subject to the same assumptions. 

1. The fluid is incompressible . 

” 2. The flow is potential. 

3. The circulation is such that, after Kutta— J oukowaki , the 

■» stagnation point occurs at the trailing edge of the airfoil. 

4. The wing is represented by a thin vortex sheet in the chord 
piano having a plan form identical to the wing plan form. 

5. All vertical displacements can be ignored. This moans, 
for instance, that (a) when camber is introduced, the chordtri.se 
angular variation is considered but not the chordwiso vortical 
displacement; (b) when angle of attack is considered no vertical 
chordwise displacements are considered; and (c) the trailing vortex 
shoot lies always in the some horizontal plane as tho wing. This 
assumption strictly limits the analysis to uncomborod wings at 
zero angle of attack; such limitations, however, can bo moderately 
exceeded. 

In replacing the wing by a vortex shoot, tho strength of the 
vorticity y at any point is related to the differential pros sure 
Ap at that point by 


&& «= pV> 
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The problem of obtaining tho loading, or distribution of Ap, 
over tho "wing is thus rosolvod into that of obtaining tho strength 
of vorticity y within tho plan form. The control condition which 
is enforced to obtain tho distribution of y is that no flow can 
occur through the vorticity sheet, or in other words, that tho 
downwash produced by the vorticity is proportional to tho slope of 
the sheet at any point within its limits. The determination of y 
would bo exact if its distribution were considered continuous and 
if the foregoing condition wero enforced at an infinite number of 
points. Such an exact determination is impractical; consequently, 
simplifying approximations must be introduced. Tho simplifications 
generally used are those cf (l) concentrating or restricting tho 
continuous vorticity chor&wiso and/or epanwise in order to maJco 
the determination of its distribution amenablo to mathematical 
treatment; and (2) representing the distribution of vorticity or 
of circulation by a mathematical expression, usually a serios, 
containing a finite number of unknown coefficients where an inf ini to 
number are generally required for exactness; and ( 3 ) limiting tho 
number of control points at which the condition of no flow through 
the shoot is satisfied. Tho differences in tho various mothods 
developed for predicting the distribution of vorticity arise, 
therefore, from (l) the manner of concentrating or restricting the 
vorticity; (2) tho differences in tho form of tho mathematical 
expressions used to describe tho vorticity distributions; and 
( 3 ) the choice in number and location of tho control points and tho 
precise mathematical procedure usod to obtain a solution. 


The Falknor Method 


Tho wing is first ccnBidorod as a- continuous sheet of 
vorticity whoso strongth distribution is expressed by tho double 
sorios : 



in which q = § and 8 


cot 

L 

2 ( a o,o 

+ iTao,i + Tp ao, 2 

+ sin 

0 (ai,o 

+ TT a l,I + 

Tf®ajL,l 

+ sin 

20 ( 02,0 

+ 

+ rf 2 a. 

*f" • • • 

a m,n T* 1 

sin m 9 


cos'~ 

1 x* 

^72 
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Evaluation of tho unknowns is performed in the follow- 

ing manner by 

1. Concentrating the vorticity both chordwise and spanwlso 

into a system of 84 finite horseshoe rorticos (fig. 1(a)) 

2. Expressing the circulation of these vortices in terms 

of the unknowns in equation (l) (appendix A) 

3. Summing at a number of control points on the wing tho 

downwash produced by all the vortices of the subject 
system and camputod by means of the Biot-Savart law 

4. Equating tho downwash angle thus determined to tho slope 

of the plate at those points thereby forming equations 
involving tho unknown coefficients 

5. Solving those equations simultaneously to eva3.ua.te the 

coefficients 

Substitution of those values in equation (l) gives the desired 
expression for the load distribution. 


The Mutterperl Method 


Mutterperl considered only spanwise distribution of vorticity. 

In such an approach "the chordwise distribution of vorticity is 
concentrated into the circulation of a lifting line. (See fig. 1(b).) 
The distribution of this circulation along the line is then repre- 
sented by the Fourier series 


r 


4kVc 0 sin a 


00 

&£n+i sin (2n+l) cp 

n=G 


( 2 ) 


The unknowns to be evaluated to obtain the distribution of T 
are the coefficients aan+i. The downwash produced at points on 
the wing by the lifting line and its trailing vertex system can be 
expressed in terms of these unknowns by application of the Biot-Savart 
law to this equation. (See appendix B.) Equating the expression 
for downwash angle to the elopes of the mean camber lines at these 
points produces a set of equations which contain the unknowns a. 2 n+i J 
simultaneous solution of these equations evaluates the coefficients. 
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The Weisslnger Method 

Fran extensions of the Multhop procedures, Weissinger developed 
two methods of obtaining span loading, one "based on lifting surface 
concepts, the other on lifting line. The lifting surface method, however 
amounted, to little more than a substitution of the theoretical ' 

additional chordwise loading, represented by y = constant x V cot 
for the concentrated load of the lifting-line method. According 
to Weisslnger the surface method proved to "be considerably longer, 
and gave results with an accuracy only slightly superior to those cf 
the line method. For this reason, only the la.tter is described 
heroin. 

As in the Mutterperl method, the continuous chordwise distri- 
bution of vorticity is concentrated into the circulation of a 
lifting line. (See figure 1(c).) The distribution of this circula- 
tion is then specified by 


G(<p) 



n=l M-i=l 


sin Picpn sin Picp 


( 3 ) 


The circulation F(y) ia represented nondimens ionally as G-(cp) 
in this expression and the unknowns to be evaluated are G^, the 
circulations at specified locations along the lino. The downwash 
produced at points within the plan form by the lifting lino and 
its trailing vortex system can bo expressed in terns of these 
unknowns through application of the Biot-Snvart law to equation ( 3 )* 
(See appendix C.) Equating the expressions for downwash angle so 
obtained to the slopes of the mean camber lines at these points 
results in a set of equations with unknowns Gn. Simultaneous 
solution of these equations evaluates the unknowns. 


EXPEP BiENTAL DATA 

Pressure data were taken at a tunnel speed of 9 ° miles per 
hour which corresponds to a Reynolds number of approximately 
9,000,000. Data were taken over an angle— of— attack range —3° to 
9 °. Plots of the chordwise distribution of pressure coefficient 
P =(p-P D )/q at several spanwiBe stations were drawn and integrated 
to obtain the local lift at these stations. These values of local 
lift were then plotted against angle of attack, and the resulting 
local lift-curve 1 slopes were used to obtain the curves of the spanwise 
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distribution of additional load and of additional lift coefficient 
shown herein. The maximum error in any local lift— curve slope as 
the result of scatter, etc., is estimated to be 0.002 per degree. 
Such an error would produce a variation of the distribution curves 
of about one-half to one— third the magnitude of the discrepancy 
between the theoretically computed and the experimentally obtained 
curves . 


RESULTS AND DISCUSSION 

Comparable spanwise distributions of the loading coefficient 
c 2 c r)/QLCav as calculated by the three methods and as determined 
from the experimental surveys are presented in figure 2. Similar 
presentations of local lift coefficient cj/Cl are presented in 
figure 3 . The theoretically predicted values of lift— curve slope 
dQk/do, and spanwise center of pressure position for the different 
wings are presented in table I. 

From figures 2 and 5 it Is apparent that all the theoretical 
methods tend to predict higher loadings at the center and lower 
loadings at the tip than were measured. In general, the Falkner 
method errs less in this respect than do the others. Mutterporl 
distribution representations for the swept-back and unswept wings 
are only slightly less accurate than those of Falkner. On the 
other hand, for the swept— forward wings the Muttorporl distribu- 
tions departed from the experimental distributions to the extent 
that they must be considered unusable. Weissinger distribution 
representations were equally accurate for swept-back and swept— 
forward wings. The average accuracy for this method was only 
slightly lees than that of the Falkner method. 

In rogard to the center of pressure position and lift— curve 
slope, the closest predictions in all instances were those made by 
the Falkner method. The Mutterperl method, in the range in which 
its applications may be considered usable, was also quite accurate. 
The Weissinger method gave good center— of— pressure positions In all 
instances and accurate values of lift— curve slopos In all instances 
except for the +45° swept wing. \ erccr -t'revn a 4"'.' 

~j\u+ ; i on . 

The time studies of tho calculations indicate that the Falkner 
method takes from 24 to 32 hours. Tho greater part of this time, 
l 6 to 20 hours. Is consumed in determining the values of the down- 
wash factor F for the different vortices. The major part of the 
remainder is needed for the solution of the simultaneous equations, 
which often prove to bo ill conditioned. Tho Mutterporl method 



NACA EM No . A7C31 


14 


takes from 20 to 28 hours, the greater part of the time "being 
consumed in the Simpson rule integration of the factors Fq 1 to F s f . 
The Weissinger method using m = M = 7, takes only 2l to 3 hours, 
in which there is no phase, that consumes an outstanding amount of 
time . v — .4 


i i ^ 


tn seY-r-i 


It has "been stated previously that the Bection lift-curve 
slope c i of all sections on all five wings was assumed to "be 
0.1030 per degree. The thickness variations from root to tip. 
however, indicate that variations in c 1 probably exist for each 

wing. Unfortunately, the distortions of the sections resulting 
from the manner in which the wings were constructed preclude an 
exact determination of what the variation might be for all but the 
unswept wing . For this reason, the readily applicable correction 
to theory (see appendixes) for a variation in ci a was not 
included in the computations. While this correction would 
account at least in part for the aforementioned discrepancies 
between theoretical and experimental loading distributions, it 
should not alter the relative evaluation of the three methods. 


In considering the three methods it should be noted that two 
of them, those of Weissinger and Mutterperl, have identical aero- 
dynamic approaches and differ, only in the mathematical procedure . 

It would be expected, therefore, if no compromise were made in the 
mathematical accuracy (i.e., if a large number of terms were used 
in the series)., identical results would be obtained. Further, if 
similar limitations were Impressed upon the two methods it might 
well be assumed that results of comparabio accuracy would be 
obtained. The failure of the Mutterperl method to predict 
acceptable loadings on the swept— forward wings is inexplicable on 
these grounds and, as a result, must be attributed to an inconsist- 
ency introduced in the mathematical development. An additional 
advantage of the Weissinger method is that it lends itself to the 
pretabulation of a number of constants which are applicable to the 
solution for any plan fozm. It Is because of this that the Weissinger 
method proved less time consuming than that of Mutterperl which 
cannot be handled in this manner. In general, then, it is apparent 
that the Weissinger method offers several advantages over the 
Mutterperl method, which, however, stem entirely from the mathematical 
phase of the solution. Insofar as the aerodynamic concepts are 
concerned neither method should be expected to be 3uporior. 
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The Falkncr method offer3 a definite aerodynamic advantage in 
that tho wing is represented by a lifting surface rather than a 
lifting line. From a consideration of only the sparcwise distribu- 
tion of loading* the time required to use the Falknor method appears 
excessive when the very minor improvement in accuracy is recognized. 
However* if surface loading or chordwise loading were desired* the 
method would undoubtedly shew marked superiority. The relatively 
long period of time required to obtain a solution by this method 
is in great measure a result of the large n um ber of purely mechanical 
functions inherent in the method. It can be expected that such 
processes are amenable to handling by mechanical means if sufficient 
use is to be made of the method to warrant their construction. One 
such aid of relatively simple form has been applied In other span- 
loading computations using the Falkner method and resulted In cutting 
the computing time by 30 percent with no soricue losa In accuracy. It 
consisted of constructing a large-scale contour chart of the down- 
wash field around a horseshoe vortex and using this in conjunction 
with an appropriately scaled drawing of the wing to read directly 
the downwash at the various control points. 

A further advantage of the FaUmer method over the lifting— line 
methods can be seen In the increased flexibility resulting from the 
system of finite vorticos which permit application of this method to 
a variety of plan forms beyond the scope of the other methods the 
lifting line pattern and control— point positions of which are 
fairly rigidly specified. In this regard, Falkner has successfully 
applied the method to a pterodactyl wing (reference 5) and to a 
wing with a parabolic 0.25c line (reference 6). It should be 
remembered* however* that should the plan form be of such a nature 
as to require a modification of tho vortex lattice, tho work 
involved will be considerably Increased. 


CONCLUSIONS 

From the results of the subject Investigation tho following 
conclusions have been drawn: 

1. Where an over-all study of the effects of sweep and plan 
form on span loading* lift— curve slope, etc.* Is desired and 
whore good accuracy Is desired for minim um effort* the Welssinger 
method Is most useful. 



1 6 


NACA BM No. A7C3I 


2 . Where a detailed study of a specific wing is desired and 
utmost accuracy is important even at the expense of considerable 
computing effort, the Falhner method should bo usod. 

3. The Muttorperl method offered no advantages over the other 
methods either in terms of accuracy or facility and honco is not 
recommended for uso. 
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TABLE I.- COMPARISON CEE THEORETICAL AND EXPERIMENTAL VALUES 
OF LIFT-CUEVE SLOPE AND SPANWISE CENTER CF PRESSURE 


Wing parameters 

Lif-t-ciArve elope dajJdjs, (deg -1 ) 

A | AS ; X 

Falkner 

— 

Mutterperl j We is singer 

Experiment 

-45.2 

-29.6 

.9 

31.0 

46.4 

2.99 ! 0.376 
4.45 j .405 
4.47; .542 
4.66 ! .442 
3.45; ,4l8 

0.0419 

.0573 

.0633 

. 0638 ^ 

.0509- 

O.G4005 

^0652 

.0615 

.0495 

0.0450 

.0585 

.0640 - 
•0^31 

, '0470' - 

___ - 

0.0422 

.0580 

.0660 

.0668 

Wing parameters 

, — 

Spanwlse center of pressure, Tin-n i 

A 

AR 

X 

Falhier 

Mutterperl 

Weissinger 

Experiment 

-45.2 

-29.6 

.9 

31.0 

46,4 

2.99 

4.45 
4.47 
4.66 

3.45 

0.376 

.405 

.542 

.442 

.418 

O .398 
.408 
.429 ^ 
.439 
.446 

0.385 

^352'T 

.426 

.434 

.438 

0.399 

.403 

.425 
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FIGURE I.- THE MANNER OF CONCENTRATING THE 
VORTICITY FOR THE METHODS OF FALKNER, 
MUTTERPERL AND WEISSINGER. 
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APPENDIX A.- PERTINENT HSFOEMATION AND COMPUTING- PROCEDURES 
FOE USE WITH TEE FALKNER METHOD 

Selection of the Yortex Pattern 

The relative strength of the circulation of the vortices in a 
network, as expressed in terms of tho unknowns in the sories 
equation (l), depends upon tho vortex pattern and tho torms in the 
series only, not upon wing shape. Tables of such circulations 
can bo set up for uso with any specified pattern. Falkner, on the 
basis of his applications, solocted the pattern of 84 vortices 
shown in figure 1(a) as suitable for most wings. While It is 
recognized that other patterns might produce moro accurate results 
in particular’ instances, tho advantages of -this regular pattern 
in reducing the computational work aro groat and hence it was 
used for all applications included heroin. 


Limitations of tho Sories 

Tho number of torms In equation (l) required to obtain a good 
approximation of the load distribution depends on the rapidity with 
which tho sorios converges for each application. For tho calculation 
of symmo tried loading, Felknor concluded that a minimum of throe 
chordwiso and throe spanwiso torms (nine unknowns) should ho usod 
for all swopt wings, while a minimum of two spanwiso and throo chord- 
wise torms (six unknowns) should bo usod for straight wings. 

It should bo roccgnizod that, as it is given, this series 
will not converge when attempting to approximate a surfaco loading 
whore discontinuities exist such as thoso rosulting from partial, 
span flapB. A slight modification to the scries, bowovor, will 
enable it to approximate tho loading wfaoro such a discontinuity 
occurs. Falkner has determined tho necessary modification In his 
investigation of wings with flaps and ailerons dof looted. (Soo 
roforonco 7-) 


Determination of Circulation of Network Yorticos 

Once tho vortex pattern and number of teiras in tho basic 
serios havo been established, the circulation of tho vortices as 
expressed In torms of the unknowns in equation (l) can bo determined 
by replacing tho continuous vorticity chordwiso and spanwise of 
equation (l) with tho concentrated stopped loading of tho lattice. 
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The chord-wise concentration of the load is determined by the 
condition that, at points located midway "between the loads (at 
one— quarter , one-half , and three-quarter chord), the downwash 
produced by the four chordwise loads IV be the same as would be 
produced by the continuous chordwise loading in two-dimensional 
flow, and the limitation that the sum of the isolated loads be 
equal to the integral of the continuous load. 

When only the first term in the chordwise series of equation 
(l) is considered. 


7 = ^ 


T a o,n 



(Al) 


or, since only chordwise loading is being considered, all factors 
which are not a function of the chordwise variable 9 can be put 
into a constant A where 


then 


A = 8s tan a<\/l— q 2 2 Tp a 0jn 

n=0 


(A2) 


7 


c 



(A3) 


Then it can be shown that if the flow is considered two dimensional 
the downwash angle at any point along the chord is 


iL _ 1 A 
T 2 c. 


and 


L, a * to 



cot & cos 0 dS 

C 


_ JtA7 
o 


(A4) 
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The following four equations may- therefore he obtained: 


1,A + ^2jA + ^3, A + *4, A - ^ TA 


(a) \ 


-8 T 2j a — ^3, A - ^4,A = « VA (h) • 
fc*i, A +^ 2j A“8l3,A -§T 4 , A (c) U5 

|ri,A + |r 2j A + er 3jA - 8r 4 , A = « va (d) J 

■where equation (a) equates the integral of the continuous loading 
to the sum of the circulations of the separate loads, and equations 
(b), (c), and (d) equate the downwaoh at the three pivotal points 
as produced by the continuous loading to that produced by the four 
loads of r v . A simultaneous solution of these equations gives tho 
increments of circulation of the four chordwise vortices which are 
equivalent to the continuous loading represented by term 1 of 
equation (l) . 


A similar solution when 


7 » sin 6 
' c 


where 


B = 8s tan a 


Tp U-Ijn 


gives the circulation— increment distribution which is equivalent 
to the continuous loading expressed by the second chordwise term 
of equation (l), and so forth. 
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The circulation 3 of a specific vortex may then he expressed by 


r v = r v ^ + r VjB +r v>C ••• 


(AS) 


or if 

r V,A " = BV«G v>b j r v,C = 

By = AVitG v>A + BVitGy^g + CVitGy^Q . . . (A9) 

The substitution into this equation of the values of A, B, etc . 
introduces the spanvise variable rj. 

r v,q = 8rtsV tan ct JL-Zf [Gv,A (a<>,o +^,1 + n 2 ao, 2 ...) 

+ Gv,B ( a l,o + n ai,i + ~ s a 1}S + . . .) 

+ GV,C ( a 2,0 + ’tT a 2 , 1 + t) 2 a 2,2 + ■••) 

+ . ..] (A10) 


Since the circulation of specific horseshoe vortices is now 
being considered, the circulation no longer varies continuously 
along the span but remains constant throughout the length of the 
bounded lines . This is equivalent to the assumption that the 
continuous loading is stepped at intervals equal to the length 
of the bounded lines of the network vertices. The continuous 
variable rf of equation (l) or equation (A10) is therefore replaced 
in a new equation by specific values n of if which indicate the 
midpoints of these lines. His new equation ■which expresses the 
circulation of any network vortex is then 


3 

The values of Gv,a Gv,B and Gv,C presented by Falkner in 
reference 1 were found to be in error. Under the direction 
of Mr. Arthur Jones these values were recomputed at Ames, 
and the values so obtained are presented in table Al. 
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r v ,u - 8rtsV tajl [g V jA (a 0j0 + H a G#1 + P 2 a c#2 + ...) 

+ G V,B ( a i,o + ^ a i,i + ^ a i, 2 + *■•5 

+ G V,C ( a 2j O + V- a 2,l + ^ a 2,2 + 

+ ... ] (All) 

or for a symmetrically loaded wing 

r v,P ** SitsT tan a'/l-u 2 (Gy^A (ao j0 + P 2 a ©, 2 + M* &o j4t + ...) 

+ G V,B K,o + & a l,s + ^ a l,4 + 

+ (^ 2,0 + P ^■s j z + P a 2 J 4- **" •••) 

+ ...3 . (A 12 ) 


Examination of this equation will show that, as has previously been 
indicated, the known parts of the equation p, Gy 3 a> Gv,Bj etc., 
are indep endent of wing shape. The products of these factors 
p n a/ 1— p 2 Gv,Aj p n a/1— P 2 Gv,Bj etc., have been tabulated for use in 
any application in which the 81j— • vortex pattern iB used. (See table All.) 


Selection of Control Points 

Since one equation iB formed at each point and since there 
should be the same number of equations as there are unknowns, the 
total number of points selected is determined by the total number of 
unknowns retained in the series equation. Further, the spanvise and 
chordwise distribution of control points must correspond to the 
number of spanwlse and chordwise terms retained in the series. The 
locations of the pointB chordwise and spanwlse are limited to posi- 
tions midway between or on the center line of the vortices. Aside 
from these limitations, the exact choice of location remains a matter 
of experience. Falkner found that for a calculation of symmetrical 
loading the arrangement presented in figure 1(a) resulted in good 
accuracy for wings with sweep. This arrangement has been used in 
all the calculations presented herein. 
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Determination of the Downwaeh 

The downwaeh produced by a simple horseshoe vortex of circulai- 
tion T is expressed by (reference 8) 


w _ 1 

V ” kit V 


F 


(A13) 


where 


F = 


_ x*-j(x*) s + (y*+l) a x*-V(x *) 2 + (y*-l) 2 ‘ 


x* (y*+l) 


x* (}■*-!) 


(Aik) 


The downwash produced by a network vortex is then from equations (Ml), 
(AI 3 ), and (Aik), 

^ = kO tan a */l-|u 2 C Gy, A ( a o, o + h &o, 1 • • • ) 

+ G v ^ b (& X} q + n a. x> x + . . . ) 

+ Gv,C ( a 2> O + a 2, 1+* • • ) 

+ ... ] F (A15) 


or for a symmetrical wing 


^ = kO tan a\Al_"p2 [ G VjA (a 0#0 + p 2 a 0 ^ 2 + . . . ) 

+ ( a i, o + a i, 2 + •••) 

+ G VjC (a 2j0 + p 2 a 2j2 + •••) 

+ . . . 1 F (Al6) 

The coordinates x*, y*, and consequently the factor F can be 
determined readily from wing geometry. In this regard, plots of the 
f unction F versus x* from 0 to 20 have been prepared at values 
of y* = 0, 2, k, 6 ... kO; however because of their size those 
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charts as such have not been included in the report, but the tabular 
data necessary for their construction are given in table AIH. In 
addition, exam i nation of (Alb) will reveal that if j* is constant 

F (— x* ) = Fi + F 2 ■- F (x*) (A17) 

where Fi + F 2 is a function of y* only 

*1 + Fa = ~r - -~r (A18) 

y*+l y*— 1 

For this reason table AIII contains only positive values of x*, 
and the function Fi + Fa is presented in table AI7. 

Summation of the downwash at any control point now ra suits in 
an expression containing the unknowns aj^ n and. their numeri cal 
coefficients which are products of the tabulated values p n ./l— y 2 , 

G v Aj G V b, G-y q, and F. In this regard it should be noted that in 
th4 summti one 3 for a symmetrical (about root chord) wing, the down- 
wash factors F for symmetrically located vortices may be added 
together prior to the multiplication of these factors by the circula- 
tions of the vortices, since in this instance the circulation of such 
a pair of vortices will be identical. 


Solution for Additional Loading 


To obtain the additional loading, the wing is considered a flat 
plate the slope at which any point is tan a = Substitution of 
thi s' value into the downwash expressions, as evaluated at the 
several control points, results in a set of equations with unknowns 
a m n* Simultaneous solution of these equations evaluates the • 
unknowns which can then be ' introduced into aquation (l) to produce 
the desired expression for - additional loading. 


The following expressions can now be derived readily from 
equation (l) : 


da 


AR it 2 
16 


(l6so o+8ai o^^o g+Ss-i 2 ^ s o 4+ a i 4 ) 


(A19> 


2l = (^) 16 yl-Tp [ 2a 0i0 -fa 1 ,o+Ti g (2a 0tg 4a :ug )+q 4 (2a 0< , 4; -fra lt ^)3 


CL 


jt [l-(l-k)Ti 1 ( l6 a 0 ,o +3a i J o +fe o,2 4£a n2^ & o,4 +a i/4 


) 


(A20) 



2k 
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o l cy = 32jl-rf [2ao .O+Q-a.o+ri lg (2a 0 . a 4-8 - 1 g^~ 4 I gga f 4 +a -1^.4.)I ( A 21) 
Cl cav Tt(l6ao,o^ai,o^^,2 +2 s-i,2 +2a o,4 +a i,4) 

t, = 3S [35(gao.o+aa.o)+l 1 )-( 2 ao. g +ax i . g )+8( 2ao J . A +a lj .4)] i (a22) 

Cp 105rt(l6a o ^ o +8a 1 ^ o +ite o ^ 2 + > 

into which idle coefficients a m n , must be substituted to obtain 
the quantities indicated. 


Solution for Basic Loading 

The determination of the basic loading on a wing with camber 
and twist can be accomplished in several way3. The simplest of these 
is to calculate tho total loading, basic plus additional, at same 
finite lift coefficient and then to subtract frcm this the additional 
loading as calculated for that lift coefficient. A solution for the 
total loading on a cambered and twisted wing is identical with that 
of a flat-plate wing up to the formation of the simultaneous equations. 
For the flat— plate wing all local geometric angles of attack wero 
identical to the wing geometric angle of attack; in this instance 
local geometrical angles of attack are in addition a function of 
the camber and twist. 

If the midwing section of the wing is choson as a reference end 
set arbitrarily at seme angle Og, then the local geometrical 
angles of attack at the various control points are known exactly; 
however, the angle a of tho reference from the zoro lift angle of 
the wing is not known. To obtain a solution undor these circumstances 
the values tan aiocal tan Oa substituted for w/V and 
tana, respectively, in the dawnwash expressions , and a solution far 
the coofficients o^n is obtained in which, however, those coeffi- 
cients will be in error by the factor tan a/tan a s . If these coeffi- 
cients and tho factor tan a are then introducod into tho expression 
for lift coefficient. 


° L = tan a ^ l6a o,c5 f ^ a i,c3 t ^ a o,^ 2a i,^ 2a o,4 +a l,4^ (^3) 

the lift coefficient for tho angle of attack a will bo obtainod, 
since tho error introduced by using ag will bo nogatod by tho orror 
in tho coefficients In other, words, the rosult is the samo 

as if the correct valuos of am,n and tan a had boon inserted into 
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equation (A23). Similarly introduction of tan as and the incorrect 
values of a m ^ n into the following: 

c l c t] ■ tanayi-rj 2 [2a 0 ^ 0 ^ 1 ^ 0 +TpC2a 0 ^ 2 +a lj2 )+q 4 (2a 0 ^+a 1 ^ 4 ) ] 

(A24) 

will result in the values of the ordinates of the total loading 
curve for tan a Now if a solution is effected for the additional 
loading., as previously described, and the value of the lift— curve slope 
dQk/da thus obtained from equation (A19) is divided into the valiie 
of Cl obtained from equation (A 23 ) , the correct value of the wing 
angle of attack tan a will result. If this value and the coeffi- 
cients n of the additional loading are then substituted into 
expression (A2^), the ordinates of the additional loading curve will 
be obtained. Subtraction of these from the ordinates of the total 
loading curve will result in the ordinates of the desired basic 
loading curve . 


Correction for Section Lift— Curve Slope 

Through the general development of the method all section lift- 
curve slopes were assumed to be the theoretical 2 k per radians 
(O.IO 96 per deg). As this assumption is not valid for all sections 
the final expression for vorticity will be in error unless a correc- 
tion is applied. If the section lift— curve slope is the same at all 
sections of the wing, the error may be corrected by simply multiply- 
ing each coefficient a^ n by the ratio of actual section lift- 
curve slope to theoretical section lift— curve slope. A varying 
section lift-curve slope can be accommodated almost as easily; 
however. In this Instance the correction must be introduced into 
equation (l) as a function of the spanwise variable q. 


Computing Instructions 

The following instructions apply to unyawed straight tapered 
swept wings without camber or twist. 

The coordinates x* ? and y*‘ relating all vortices and 
control points to the center section leading edge of the wing are 
calculated on form A(l) using the relations 
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y* f = 20 p 

= 20 |m! tan A— [1 - (l-X.) |p|] (A25) 

for the vortices, and. 


y*‘ = 20 t] 

x *' = 20 l^i "ten A - [1 - (I-*-) |nj] (A 26 ) 

for the control points. Since the wing is symmetrical x*' will "be 
the same for similarly located vortices on each wing half, and the 
values of y* f for the left wing will he the same as those for the 
rigjht wing, although of opposite sign. For this reason these values 
only need he computed for positive values of rf or p. 

The values of x* and y* relating a control point to each of 
the vortices are obtained hy subtracting the values of x* 1 and 
y* 1 of the vortices from thoee of the control point, column 9 or 3, 
form A(l) from a value in column 18 or 12, form A(l), respectively; 
x* and y* are then tabulated on a form A (2), using a separate 
form for each control point. It should he noted that since the 
coordinates of the vortices at p = 0.9625 are based upon a unit 
length y v one-quarter normal size, x* and y* for these 
vortices are four times the normally calculated values. Lastly, hy 
virtue of symmetry of plan form, the coordinates can he tabulated 
so that two value s y^* and y R * exist for every value of x*. 

These coordinates are now used to enter charts of the downwash 
function F as prepared from the values in table AIII. The values 
obtained for the vortices at p = O .9625 should.be multiplied by 
four. Because of symmetry of loading, Fg and F^ can be and are 
added together. - 

The simultaneous equations set up in tabular form in form A ( 3 ) 
are now obtained as follows: Considering the first equation or 

column 1, the second number, the numerical coefficient of ao.Oj 
is obtained by multiplying the values of Fg + F^ in column 7, 
form A(2) as determined for control point 1 by the values listed 
under a 0 0 in tablo All and summing the products. Similarly, 
the third '’number in column 1, the numerical coefficient of 0 , 
is obtained by multiplying the values in column 7 "by the values 
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listed under a a 0 in table All and summing the product. The procesB 
is repeated using the values listed under a 2 0 , a Q a , etc., until 
the entire equation is obtained. J ~ 

The second equation, column 2, form A(3), is set up in the same 
manner except that the values in form A(2) as determined for a 
second control point are used. The procedure is then repeated until 
the nine required equations are formed. 

The constant numbers, row 1 of form A(3)> are obtained as 
follows: The dovmwash at control point 1 is 

5L = kO tan a X column 1, form A(3) 
or 

(-i) = column 1, form A(3) 

V tan a,y 

Equating ^ to the slope of the plate, tan a, 

—■ = 0.0250 = column 1, form A(3) 


In like manner the constants for the other equations are also 0.0250. 

The equations are set up in this manner to facilitate their 
solution by the method outlined in reference 9* Of the various 
methods for solving a large number of simultaneous equations which 
were tried, the method of reference 9 was found to be most rapid and 
straightforward where only standard computing machines were available. 
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TABLE AI.- OEQRDWISE FACTORS G FOR VORTEX 
PATTERN UTILIZING FOUR CKGRDWISE VORTICES 


V 

%,A 

g v,b 

G v,c 

0.125 

» 

0.27337 

C. 04 902 

0.C7282 

.375 

.11680 

. 07598 

.03823 

.625 

.06947 

.07538 

-.03823 

.875 

.04036 

.04902 

. - . 072 82 
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TABLE All.- PRODUCT OF OHOHDWISE A HD SFAHWISE FACTORS FOR t m , a TO BB 
USED WITH TEE 84-VOETEX PATTERS' SYVMETRIOAL LOAD ISO 


NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 


X OVJL 
for 
*0,0 

JT- v? 

X Ov,B 
for 
*1,0 

for 

* 2.0 

nVl-w* 
x OVJL. 

for 

*o,a 

n^A-n a 
* Gv.B 
for 

*i,a 

X Ov,C 
fer 
*a,a 

fdr 

* 0,4 

x Gv,B 
for 

X 0v,C 
for 

«a,4 

0.2735? 

0.04902 

0.07282 

0 

0 

0 

0 

0 

0 

.27200 

.04877 

.07246 

.00271 

.00049 

.00072 

.00005 

0 

.00601 

.26785 

.04805 

.07155 

.01072 

.00192 

.00285 

.00044 

.00008 

.00012 

.26077 

.04676 

.069*7 

.02548 

.00421 

.00626 

.00210 

.00056 

.00066 

.25054 

.04493 

.06674 

.04008 

.00719 


.00642 

.00115 

.00171 

.25674 

.04245 

.06506 

.05918 

.01061 

.01677 

.01479 

.00256 

.00594 

.21870 

.05922 

.05826 

.07873 

.01412 

.02097 

.02855 

.00508 

.00766 

.19521 


.05200 

.09665 

.01716 

.02548 

.04688 

.00841 

.01249 

.16402 

.02941 

.04569 

.10497 

.01882 

.02796 

.06719 

.01206 

.01790 

.11916 

.02157 

.05174 

.09655 

.01731 

.02671 

.07818 

•01402 

.02085 

.07417 

.01550 

.01976 

.06870 

.01252 

.01850 

.06564 

.01141 

.01695 

.11680 

>07598 

.0582$ 

0 

0 

0 

0 

0 

0 

.11622 

.07560 

.05804 

.00116 

.00075 

.00058 

.00001 

.00001 

0 

.11444 

.07445 

.05746 

.00458 

.00298 

.00160 

.00019 

.00012 

.00006 

.11142 

.07248 

.03647 

.01005 

.00665 

.00528 

.00090 

.00069 

.00029 

.10706 

.06964 

.03504 

.01712 

.01114 

.00660 

.00274 

.00179 

.00090 

.10115 

.06580 

.05311 

.02529 

.01645 

.00828 

.00652 

.00411 

.00207 

.09544 

.06078 

.05058 

.03364 

.02188 

.01101 

.01211 

.00788 

.00396 

.08541 

.05426 

.02730 

.04087 

.02659 

.01358 

.02005 

.01305 

.00666 

.07008 

.04669 

.02294 

.04486 

.02918 

.01468 

.02871 

.01868 

.00940 

.06091 

.05512 

.01666 

.04124 

.02683 

.01350 

.05340 

.02173 

.01095 

.05169 

.02061 

.01057 

.02954 

.01909 

.00961 

.02709 

.01769 

.00890 

.06947 

.07698 

-.05825 

0 

0 

0 

0 

0 

0 

.06912 

.07560 

-.05804 

•00069 

.00075 

-.00038 

.00001 

.00001 

0 


.07445 

-.03746 

.00272 

.00298 

-.00160 

.00011 

.00012 

-.00006 

.06627 

.07248 

-.05647 

.00597 

.00665 

-.00328 

. 00055 

.00069 

-.00029 

.06567 

.06964 

-.05504 

.01018 

.01114 

-.00660 

.00165 

.00179 

-.00090 

.06016 

.06580 

-.03311 

.01504 

.01645 

-.00828 

.00376 

.00411 

-.00207 

.06558 

.06078 

-.05068 

.02001 



.00720 

.00788 

—00596 

.04961 

.06426 

-.02730 

.02451 

.02659 


.01X91 

.01305 

-.00656 

.04168 

.04559 

-.02294 

.02668 

.02918 

H3E im 

.01708 

■2 i.m 

-.00940 

.05028 

.03512 

-.01666 

.02453 

.02685 

-.01560 

.01987 

IMS 

-.01095 

.01885 

.02061 , 

-.01037 

.01746 

.01909 

-.00961 

.01617 


-.00890 

.04056 

.04902 

-.07282 

0 

0 

0 

0 

0 

0 

.04016 

.04877 

-.07246 

.00040 

.00044 

-.00072 

0 

0 

-.00001 

.05954 

.04805 

-.07155 

.00158 

.00192 

-.00286 

.00006 

.00008 

-.00012 

.05850 

.04676 

-.06947 

.00547 

.00421 

-.00626 

.00051 

.00058 

-.00066 

.05699 

.04495 

-.06674 

.00592 

.00719 

-.01068 

.00095 

.00115 

-.00171 

.05495 

.04245 

-.06506 

.00874 

.01061 

-.01677 

.00218 

.00265 

-.00394 

.05229 

.05922 

-.05826 

.01162 

.01412 

-.02097 

.00419 

.00508 

-.00755 

.02882 

.05501 

-.05200 

.01412 

.01715 

-.02548 

.00692 

.00841 

-.01249 

.02422 

.02941 

-.04569 

.01550 

.01882 

-.02796 

.00992 

.01206 

-.01790 

.01759 

.02157 

-.03174 

.01425 

.01751 

-.02571 

.01164 

.01402 

-.02083 

.01095 

.01530 

-.01976 

.01014 

.01252 

-.01850 

.00940 

.01141 

—01695 










TABLE AIII DOWMWASH FACTOR F IN THE FIELD OF A HORSESHOE VORTEX 

[At positive values of x* only] 


% 

F + (2/x) 

F 

f" = 0 

y* = 2 

y* = 4 

6 

8 

10 

12 

14 

16 

18 

20 

n 

O 

2.00000 

—.66667 

-.13333 

-.05714 

-.03175 

-.02020 

-.01399 

-.01025 

-.00785 

-.00619 

-.00502 

.1 

1.90020 

-.62234 

-.12976 

-.05616 

-.03135 

-.OI999 

-.01387 

-.01018 

-.00780 

-.00616 

-.00499 

.2 

1.80200 

-.57875 

-.12625 

-.05519 

-.03093 

-.01980 

-.01375 

-.01011 

-.00774 

-.00612 

-.00497 

.3 

1.70647 

-.53652 

-.12271 

-.05421 

-.03053 

-.01959 

-.01364 

-.01003 

-.00770 

-.00609 

-.00493 

.4 

1.61484 

-.49621 

-.11920 

-.05326 

-.03011 

-.01938 

-.01352 

-.00996 

-.00765 

-.00605 

-.00491 


1.52788 

-.45819 

-.11573 

-.05226 

-.02973 

-.01919 

-.01340 

-.00988 

-.00760 

-.00602 

-.00489 

.6 

1.44603 

-.42269 

-.11228 

-.05130 

-.02934 

-.01897 

-.01329 

-.00981 

-.00755 

-.00599 

-.00486 

.7 

1.36954 

-.38981 

-.10889 

-.05034 

-.02893 

-.OI878 

-.01317 

-.00974 

-.00750 

-.00595 

-.00484 

.8 

% 

CVJ 

1 — 5 

-.35957 

-.10955 

-.04938 

-.02854 

-.01857 

-.01305 

-.00967 

-.00745 

-.00592 

-.00481 

.9 

1.23254 

-.33186 

-.10227 

-.04843 

-.02814 

-.01838 

-.01293 

-.00959 

-.00740 

-.00589 

-.00478 

1.0 

1.17158 

-.30655 

-.09904 

-.04750 

-.02774 

-.01817 

-.01282 

-.00952 

-.00736 

-.00584 

-.00476 

1.5 

.92963 

-.21017 

-,o 84 oo 

-.04292 

-.02581 

-.01718 

-.01224 

-.00915 

-.00711 

-.00568 

-.00464 

2.0 

.76393 

-.14956 

-.07093 

-.03863 

-.02394 

-.01620 

-.01167 

-.00880 

-.00686 

-.00551 

-.00451 

2.5 

.64594 

-.11032 

-.05986 

-.03468 

-.02214 

-.01526 

-.01112 

-.00844 

-.00663 

-.00534 

-.00439 

3.0 

.55849 

-.08398 

-.05066 

-.03107 

-.02045 

-.01435 

-.01057 

-.00810 

-.00639 

-.00517 

-.00427 

4 

.43845 

-.05255 

-.03683 

-.02492 

-.01739 

-.01264 

-.00953 

-.00743 

-.00593 

-.00485 

-.00403 

6 

.30574 

-.02555 

-.02100 

-.01632 

-.01255 

-.00974 

-.00770 

-.00620 

-.00508 

-.00423 

-.00357 

8 

.23444 

-.01489 

-.01318 

-.01112 

-.00918 

-.00752 

-.00620 

-.00515 

-.00433 

-.00367 

-.00314 

10 

.19012 

-.00969 

-,00893 

-.00791 

-.00685 

-.00587 

-.00500 

-.00428 

-.00367 

-.00318 

-.00277 

is 

,12889 

-.00438 

-.00422 

-.00397 

-.00368 

-.00335 

-.00305 

-.00274 

-.00247 

-.00223 

-.00200 

20 

.09750 

-.00248 

-.00243 

-.00234 

-.00223 

-.00211 

-.00199 

-.00185 

-.00171 

-.00157 

-.00148 


u> 

ro 
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TABLE AIII.- Concluded 


f 


■ 

r 

B 



y * 22 

24 

26 

28 

30 

32 

34 

36 

— 

38 

40 

0 

-.00413 

-.00348 

-*♦00296 

-.00256 

-.00222 

-.00196 

-.00173 

-.00154 

-.00139 

-.00125 

•i 

-.00412 

-.00346 

-.00295 

-.00255 

-.00221 

-.00195 

-.00172 

-.00154 

-.00139 

-.00125 

.2 

-.00410 

-.00345 

-.00294 

-.00254 

-.00220 

-.00195 

-.00172 

-.00154 

-.00138 

-.00124 

• 3 

-.00409 

**• 00343 

-.00293 

-.00253 

-.00219 

-.00194 

-.00171 

-.00153 

-.00138 

-.00124 

.4 

-.00407 

-.00342 

-.00291 

-.00252 

-.00218 

-.00193 

-.00171 

-.00153 

-.00137 

-.00123 

.5 

-.00405 

-.00340 

-.00290 

-.00251 

-.00218 

-.00192 

-.00170 

-.00152 

-.00137 

-.00123 

• 6 

-.00403 

-.00338 

-.00289 

-.00250 

-.00217 

-.00192 

-.00170 

-.00152 

-.00137 

-.00123 

►7 

-.00401 

-.00337 

-.00288 

-.00249 

-.00217 

-.00191 

-.00170 

-.00151 

-.00136 

-.00122 

.8 

-.00400 

-.00336 

-.00287 

-.00248 

-.00216 

-.00190 

-.00169 

-.00151 

-.00136 

-.00122 

*9 

-.00398 

-.00334 

-.00286 

-.00247 

-.00215 

-.00190 

— » 00168 

-.00150 

-.00136 

-.00122 

1*0 

-.00396 

-.00333 

-.00285 

-.00246 

-.00215 

-.00190 

-.00168 

-.00150 

-.00135 

-.00122 

1.5 

-.00387 

-.00325 

-.00279 

-.00242 

-.00211 

-,00186 

-.00166 

-.00148 

-.00134 

-.00120 

2.0 

-.00377 

-.00318 

-.00273 

-.00237 

-.00208 

-.00183 

-.00163 

-.00146 

-..00132 

-.00118 

2.5 

-.00368 

-.00311 

-.00268 

-.00233 

-.00204 

-.00180 

-.00161 

-.00144 

-.00130 

-.00117 

3.0 

**• 00358 

00304 

— . 002 63 

-.00228 

-.00200 

-.00177 

-.00158 

-.00142 

-.00128 

-.00115 

4 

-.00339 

-.00291 

-.00251 

-.00219 

-.00193 

-.00172 

- i 00153 

-.00137 

-.00124 

-.00113 

6 

-.00305 

-.00263 

-.00230 

-.00201 

-.00179 

-.00160 

-.00143 

-.00129 

-.00117 

-.00106 

8 

-.00273 

-.00237 

00209 

-.00186 

-.00165 

-.00148 

-.00133 

-.00121 

-.00110 

-.00101 

10 

-.00242 

-.00214 

-.00190 

-.00170 

-.00152 

-.00137 

-.00124 

-.00113 

-.00103 

-.00096 

16 

-.00180 

-.00164 

-.00148 

00 X 34 

-.00123 

-.00113 

-.00103 

-.00095 

-.00088 

-.00081 

20 

-.00135 

-.00126 

-.00116 

-.00106 

-.00099 

-.00092 

-.00086 

-.00079 

-.00074 

-.00069 


u> 

U) 
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TABLE A IV.- AUXILIARY FJNCTICN F, + F 2 FOR DETERMINING 
DOWNWASH FUNCTION F AT NEGATIVE VALUES OF x* 
[F(-x^ = F 1+ F - F(+x*) ] 


y* 

P l+ P a 

y* 

P l+ P a 

0 

4.0C00 

24 

-0. 0070 

2 

-1.2333 

26 

-.0059 

4 

-.2667 

28 

-.0051 

6 

-.1143 

30 

-.0044 

8 

-.0635 

32 

-.0039 

10 

-.0404 

34 

0035 

12 

-.0280 

13 

-.02381 

14 

-.02 05 

37 

-.00292 

16 

-.0157 

61 

-.00108 

18 

-.0124 

93 

-.00046 

20 

-.0100 

117 

-.00029 

22 

-.0083 

141 

-.00020 
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APPENDIX B — PERTINENT INFORMATION AND COMPUTING 
PROCEDURES FCE USE WITH THE MJTTERPERL METHOD 


A lifting line used to represent a wing Is placed in a posi- 
tion corresponding to the quarter— chord line of the wing. The 
distribution of circulation along the lifting line is expressed by- 
equation (2). 

No generalization can be made as to the number of terms which 
must be retained In the series to ensure acceptable accuracy. 
Mutterperl implies that four are sufficient and utilizes this number 
in all applications. It should be noted that all loadings predicted 
by the series as it stands will be symmetrical. In addition, 
equation (2) cannot satisfactorily approximate a curve containing 
discontinuities such as would be produced by flaps or ailerons. 
Mutterperl made no comment as to additions or alterations to the 
series which would enable circumvention of these limitations. As 
a result, while it is believed that such modifications could be 
included, it Is not known to what extent they would increase the 
complexity of the mathematical evaluation. 

Since it can bo shown that in a theoretical approach using a 
lifting line at the quarter— chord line, the downwash angle at the 
three— quarter— chord line most closfly approximates "the true angle 
of attack of the wing, the control points were placed along this 
line. The number of points required is dictated by the number of 
unknown coefficients retained in equation ( 3 )* The location of 
these points spanwiso *ie apparently arbitrary; however, since 
Mutterperl placed them on the right wing half at q = 0.174, O. 5 OO, 
O. 766 , 0.940 = 80°, 60°, 40°, 20°), this arrangement has been 

followed in all applications presented herein. 


Determination, of Downwash 

The expression for downwash at a control point, as determined 
through the Biot-Savart relation, is 
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l C2°+Ds=n-i -j&tm ~ Bl ,A+S%) 

n=0 

it sin (2n*l)>ir + s r cos (2n+l)cp cos cpdfj? 
s r cos A sin iff B^Jq F 

- i(i ■ * ^)f A ♦ *0- * ^ ■ »- *) 


/ . \ pnfscoB (2n+l)cp( coe <p — J>± f — j cup *1 

Zaa+i -I* 1 1/ — g1 -^ 

^ Mo F(C+E) J 

L. f^ -ssa. ffi c°g fan-i.) opdtp, + _ tan a) f^ Z po^JgniiUggffi 

2 LJn & BlVBl /Jo g 


x 


b 2 l 


+ L_ 

s 


_ A _ (At.— y* ) tan A \ p^ s cos (2nH-l)cpd<p ^ ( B1 ) 

r \ B L JJ 0 cos <? + g -■ 


Equation (Bl) may "be reduced to a simple expression, centain— 
ing only the unknowns a 2 n+x and their numerical coefficients 
simply by the introduction of wing geometry and the geometry of a 
control point. Such a redaction should therefore be carried out 
at each of the -points. 


Solution for Additional Loading 

Since for such a solution the wing is considered a flat plate, 
all geometric angles of attack became sin cc. and the constant 
factors in equation (Bl) become one. Simultaneous solution of 
these equations then evaluates the unknowns a 2 n+i, which are 
subsequently introduced into equation (2) to produce the desired 
expression for additional loading.- The unknowns can also be intro- 
duced into the following expressions derived from equation (2), to 
obtain the values indicated. 



MCA EM Wo. A7C31 


37 


dC T , 2it2bC n _ kr.s&! 
do, “ S Sa = 1+X 


(B2) 


ci _ 2(l+X) 

Cl icaj.fl— (l-X) cos cp] 


(ai sin cp + a 3 sin. 39+as sin pep +aer sin 7cp) 

(B3) 


£I2H — = -it_ (a^ sin cp+ag sin 39+% sin 59+a^ sin 7<p) (B4) 

Cl c a v j caj. ^ 


■Hep 


k_ 

3* 



3 s * _ g?_ + 
5ai 7&i 



(B5) 


Solution for Basic Loading 

As In Falkner, the basic loading on a twisted and/or cambered 
airfoil can best be calculated by determining the total and addi- 
tional loading at seme finite lift coefficient and subtracting the 
latter from the former. The procedures Involved are parallel to 
those of the FalTmer method as well. An arbitrary angle of attack 
<Xg can be selected for the root section of the wing, from which 
all local angles of attack can be measured. Substitution of 
Bin cciocal for w/V and sin a s for sin a in the expression (Bl) 
will result In a set of equations which may be solved simultaneously 
for the values of the coefficients agn+i- If aj. and sin Og are 
are then introduced into the following, the correct value of the 
lift coefficient for . the wing at the described attitude a will 
result 


C L = 4itgai sin a 
1+X 


(B 6) 


In addition. If the values of a^+x and sin a 3 are substituted 
into the following, an equation for the erdinates of the curve of 
total loading on the wing at a results . 


cic-q = (ax sin cp+a 3 sin 39+&s sin 5^+ar sin 7<p) 


V" i / 
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IT the values of coefficients obtained for a solution for the addi- 
tional loading are then substituted into this expression, the addi- 
tional loading at this coefficient ■will result. Subtraction of the 
ordinates ctf this last from these of the total loading curve will 
give those of the desired basic loading curve . 


Correction for ci a 

As in tho mothod of Ballmer, the error introduced into the 
solution by the assumption that all section lift-curve slopes c i 
were 2 it can be readily eliminated. I? the actual ci does 
not vary across the wing, the coefficients should be multi- 

plied by the ratio of actual cj a to Sir . If the actual value does 
vary along the span, this ratio Bhould be included in equation (3) 
as a function of the spanwise variable <p. 


Computing instructions 

These instructions apply only to unyawed wingB devoid of camber 
and/or twist. 

If the local angle of attack sin a is introduced into 
expression (Bl) in place of the downwash ratio w/V, this expres- 
sion can bo written 

3 

1 = ^ ( 2n+l ) a 2n+ i (F x +F 2 +F 3 F 3 f -F 4 F 4 1 +F 5 F 5 l +F 6 F 6 1 +F 7 F T «+F a F 8 1 } (b8) 

n=0 

where 

Fl = A B A L. == N ) 

2n+l BLa/A^L+BSl/ 


F e = Z£ sip (20+1)#, (B 9 ) 

s 1 cos A sin ijr 


F = 

£ 3 
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Computing form. B(l) is used to calculate all factors which remain 
constant throughout the sumnations for any one control point: Ag, 

Bg, Al, Bl, s', C, F s , F 4 , F 5 , F 6 , F 7 , F s . The computation form 
for those factors, which vary with n only throughout the summations 
for any one control point, is presented in form B(2) . Computing form 
B(3) is used to apply Simpson's rule to the integration of F 3 1 , R ; 1 , 
F 5 l , F 6 ', F 7 ‘ and F e 1 . The factor Kgn+d, which is independent of. 
wing shape and so can he applied to all wings, was calculated from the 
relation 


Ksn+i 


/-eir 

1 


cos( 2 n+l) dg> 
cos <p — cos 


(BIO) 


and is presented in the following table: 


Control 

point 

n 

1 

Ksn+i 

1 

0 

2.54928 


1 

8.26392 


2 

8.87752 


3 

6.10976 

2 

0 

2.47976 


1 

4 . 5S64Q 


2 

- 1.90944 



- 4.63656 

3 

0 

2.33120 I 


1 

0.47920 


2 

-3.47696 ; 


3 

3.39776 I 

4 

0 

2.00112 


1 

-1.92340 


2 

1.38068 


3 

-0.53112 

! 


The results of the integration are presented in form 3(3) . 
Form B(4) is the form in which the components F are multiplied 
together and the results are summed producing four equations, 
one for each control point. Lastly, the form for the simultaneous 
solution of these equations (reference 9 ) and the resulting values 
of the unknowns are presented in form B(5) ■ 
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A_= ~ /% tany\_= - 1,0*16 it r\J\.= -.22l£2 co&7\_= . 90 4 (.3 umA.* ~ 12221 cosZA. =-.£2±lS 


k’ Oli (/“A>- '±±1 M= %j31— s'- ^ioi/L = h. 2 co ^ 


i 
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9 
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1 1 
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18 
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a 
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|-<5) 

2 

<T 
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(•<1 
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► 
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20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

3 ( 

32 

33 

34 

35 

36 

68) 

CIO) 

c 

s' 

(I2) 1 

r 3 _ 

do) 

Czi) 

+toiA 

fe) 

do) 

Ft 

fcsT 

** 

C3>s # 

6) 

-(zs) 

i+ . 
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(IO) 

0^5' 

(zt) 

feft) 

F s 

s' 

(IS) 1 - 

n 

(4) 
C IZ) 

C3l) 

-CunA 

iiii 

(IZ) 

f 7 

<H) 

-Cis) 

i- 

NDuiA 

(jz) 

is 

(,W9*S 

j£gSg 

t.Wft 

74/6*6 

tt.trs?* 

Aoaa/oSi 



t.o/Qct 
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Form B‘~ Computing form 

(Underscored numbers 
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COMMITTEE FOR AERONAUTICS 

FOR MUTTERPERL’S METHOD 
are sample calculations J 
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Form B :_ Continued 
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APPENDIX C DEVELOPMENT AND COMPUTING PR0CEDUBE3 
FOE USE WHS THE WEISSIEGER METHOD 


In tihe German reports available, which described this method, 
the mathematical development was not complete . It was thought 
advisable therefore to Include the development in the present paper. 

From lifting-line theory the equation of downwash at the one- 
quarter chord line of a straight wing' is given by 



= x W g rtqp 
k-xj. b/e y-y 


<*y 


(Cl) 


This integral equation is solved by Multhopp by an Integration 
formula (reference 10) . 

The equation of downwash at any point xy of a straight wing 
is given by 

Vx y * -i- r 1 + 1 ■ r * (y) d? (C2) 

^ J-b/£ y-y { x J 


Weissenger divides this Integral into two integrals, one of which 
is the same as equation (Cl) which he solves by Multhopp* s method; 
and the other which he solves by a method analagous to that of 
Multhopp . 

The mathematical development Isas follows: 


With G = f/b V, q = 2y'/b , If = 2y/b,. ar = b/c^ 


and setting x equal to the distance to the three-quarter chord line 
x = Ctj/ 2, the equation (C2) becomes 


at „ JL T 1 • 
v 2it iJ L 1 tHT 


1 +*/l +(ar) £ (ri-^ - ) 2 G* (Tf)dW 
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or 


1 - h£l al $r 1 + Wijll W [ar(n-Ti)] 0 Yn)«T ( 03 ) 


where 


X-A=0 [ar(ri-Ti)] = - 


(C4) 


or 


La=0 (v,h) = L[ar(cos cpv - cos <^l)] = 


j " ' ' 1 p 

Vl+ar v (cos cpy — cos cpp.) — 1 
ar v (cos cpy — cos qajj. } 


The first integral of equation (C3) can be written as a function 
of cp. 



COS cp — COS (Pv 


(C5) 


whore 


COS <Py 


= COS = 2z. = T) 

m+1 b ' 


An integration formula gives 


F 1 f ^ f(ri n ) sin 

J-i m+1 L 

n=l 


Cc 6 ) 


whore cpjj = f(T)n) is the value of f(rj) at t| n . 
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Equation (C6) holds exactly if f(rf) can be represented by • 


2m— 1 


2m 

V 


fCq) = yi-n 2 ^ Avn V = 2 ft V sin ^ 
V=0 V=1 

or, in addition, 

f(Cp) = £ sin l^iCpri sil1 t 1 ! 5 ? 


n=l 


tii=l 


letting G-(cp) = f(<p) 
then 


m m 


G(cp) = ^ Gq sin UxC^ sin iiicp 


a 2 

dtp m+1 


n=l m=l 
m di 

X X 

n=l M-i=l 


Hi sin.P-x'Pn 2 ? . c0 . s 


and -the integral in equation (C5) beccmes 


(C7) 


Nov 


2 

m+1 


J it m 

G_ r V ' n. Sin [XjCfa cos mg, dqj 
Jo Z_. cos cp - cos cp v 


n=l ixi=l 


r cos ncp ^ _ it sin ncp Q 

cos cp — cos cp 0 sin cp Q 


(C8) 


or 


m . 

Z' ^ Cm-1 Sin liiCpq) 


COS UtCP 


.cos cp — COS Cpv 


dm = It V ^ s — HaShl, fcL*P.v 
J sin cp v 




Ri=l 
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so that equation (C5) becomes 



di sin UiCpn sin Picpy 
sin cpv 


lii sin sin Hitpy 

sin cpv 


(C9) 


The |ii series is independent of wing geometry and may be put into 
a coefficient usable for all wings . 


From equation (C9 ) , for d=<v 3 let 


bv,V= 


. X. 

(m+l) sin cpv 


m 

I 

Pi=l 


|i a sin 2 Picp v 


(CIO) 


and for n 4 v > lot 


m 

-n * gaffer s I “ iaIn ^ sln ^ 

dl=l 


Then equation (C5) becomes 


(cn) 


bv^vGv 



(C12) 


(Note: The summation prime indicates that the term of n = V 
should not be included.) 


Expression (C12) gives the induced angle of attack on the one— 
quater— chord line 3 at the span station v, in terms of the summation 
of n spanvise values of the dimensionless circulation. 
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Equations (CIO) and (Cll) may be simplified with the aid of 
the summation formula 


m 

V 

L 

di=i 


l-i i cos M-iX = real part o: 


<1 


m .* 111 'z 
/ 


di=l 


- m cogfm+l'x - fm+l'icos mi + 1 
2(ccs r— l) 


OT 


^v.v 


= — JB±1 

> v 4 sin ?V 


"bv. n = 


sin tPn 


(cos cp n - cos qpv) 1 




n-v 




(C13) 

(Cl*) 


The second integral of the integral equation (C3) or 




lA(V,!a)G*(Tl)dq 


(C15) 


is solved in an analogous fashion. 
The integration formula gives 


^T rt f(o)dq) - f( *)] (016) 

P-=l 


where 



HACA EM Wo . A7C31 


Dimensionless circulation is given by- 


Lett ing 


m m 


G(ep) - ^ G n V sin ^i^n sin 


n=l (.(.1=1 


G l (<p) = -2- V G n Y Hi3in Hi o n cos Hi<p 

m+1 Z_ Z_i 

n=l Hi=l 


then 


i-J?) = ^ ^i 3111 l-^'Pn 008 

Hi=l 


■£ 


G f (cp) = \ Gn f n (<p) 


Expression (C15) say be modified, so that 


±.£ X L A (v f[ x) G‘(Tf)dTi = - L A (v,'i)G l (cp)dcp 


m 
V 
= > 


: r rL A (v,u)Gn f n (<p) 
Jo L 


Then applying equation (Cl6)j ^ y La(v,|jl) G-' (rj‘)dr[ 


1 V jl (r^hAy.!?.)... 

2 it M+l t 


(C17) 


M 

V~ 

+ ) L A (v, 
H=l 


Gn — — J Vj.M+3.) ,3 (f n ,,M+i),. + ^ LA(v,p)f n p 1 
2 (M+l) L 2 U 3 J 
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2Tor simplicity let f n ^ where fj^^) is from 

equation (C17) for cp = = pit/M+1 

Then if 


6v ^ 2(M+1) 

expression (C15) becomes 


jL_ ^i L A ( v ;° 2 1 ( f n, l(£aJ£al + ^ LA(v,.u)f njl ij- 

p.*l 


(Cl8) 


m 


^ G n g v ^ n 


(C19) 


n=l 


The solution of equation (C3) is two timeB equation (C12) plus 
ary times equation (C19) or • 


m j m 

= 2 tVjVGy — ^ "bv,n GfiJ + V ar v Sv,n 


n=l 


n=l 


' el,» m t 

= 2 b v ^ v G v + ar v Sy^ v Gy + ^ ®v,n ~ ® \ ^v,n^h 


n=l 


n=l 


so that 


m 


- i 


_ V 

L.- 

n=l . ■ 


^ = (2bv,v + ary 6v,v) Gv — 2_ C^Vj.n - ^'v Sv,n) °n 


or 


f = b* v Gy 




^*Vjn Gn 


n=l 


(C20) 
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whore 


b*v ** 2by^y + ®^v 8v,v 

= 2b^ n — ar v 6v,n 


(C21) 


The prime on the summation indicates that the term of n = v 
should not he added into it. 


Equating the downwaeh to the local angle of attack of the plate, 

m 

Cty = b*y Gy — V' b*y n On (022) 


n^J 


For a. swept wing the equation of downwaeh at any chordwise 
point for y>0 is given by 


w xy 


= _L r 1 (y) 

Mb / 2 y-y 
+ 5?/b/a r(?) 


1 + 


x — jy 1 ton A 

V (x - [y | tan A) 2 + (y-y) 2 J 


Ay 


x — y tan A 


[(x - |y | tan A) 2 + (y~y) 2 ] 3 /a 


Ay 


x - y tan A j— 

[(x - jy j tan A)^+ (y-y) 2 ] 3/ ^ 


(C23) 


The first integral is the downwash duo to the trailing vortex 
sheet and the last two integrals represent the velocity induced 
by the lifting line . By integrating the last two integrals by 
parts and rearranging into dimensionless quantities, the procoding 
integrals may be put in a form similar to equation (C3) . 

‘W "f + Sr £[ Mi-to O' ft) d n (raw 


(Hote that in tho following equation (C25) the squares under the second 
radical of IiA(q,q) for rf<0 aro summed. In roforonce 3 these 
were erroneously shown as a difference . ) 



where for q<0 


t.fti.Tf) = — r* ffilt-frE- 3M^J- g (.Tlr-jl ).^ _ jl + ^BA/Il+n(ar) tan A1 g ± (or) 2 ^ | 

A ar(Tj-rf) L 1 + 2q (ar) tan A J ' ‘ 


and for q>0 


1 + 2i) (ar) tan A 


l a (t},ti) = + ar~tan a. ,(n-Jnl.U 5 t, .(ar) g C n— ri) g ~i, 

ar (tj-tj) 


}(C25) 


Equation (C24) with La(t]jtD as given by equation (C25) is for y>0 and will give 
values for only the positive span stations. For values with y <0, expressions similar to 
equations (C23) and (C25) must he derived. Values of L A (q,iy) "are needed for y<0 if an 
unsymmetrical wing is to analyzed. ' “ 

Equation (C25) with rj => cos <p v and Tf = cos and or at span station v as ar v 
for cos cp v >0 becomes, for cos cpjj. <0, 


L / v a 1 fj [ 1+ary tan A (cos ^ ~ [cob <Pu )) ] (ary) a (cos cpy - cos cpu) : 

A' 3 * ary (cos cpy “ cos $. 1 ) [ 1-2 ar v cos cp v tan A 


+ 2 tanv >r [l+ar v tan A cos q?-u3 g + (ar v ) 2 coo 5 '^ 
1+2 (ary) cos q)y tan A 


and for cos >0 


L A (v,n) = '/[1+&JV tan (cos cp v - jcoa + (ary) H cos cpy - cos q\x) k -1 

ar v (cos q)y -■ cos <tyj 


>(026) 




vn 

H 
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Equations for determining Gy at spanwise points of 

^ =» cos for any wing may therefore he determined from 

b m+1 . 

equation (C20) . 

Summarizing the computations, the relation equating dovnwa.sh 
to the local slope of the plate at m points along the span giving 
m equations with m unknowns, G^ is 

ay = b*y Gy — ^ b^y^n Gjy ^ = 1*2, . . . m 
n=l 


where 


ay = angle of attack at span station V 


— t= cos 
b/2 


m+l 


b*y = 2by^y + aTy gy^y 
b *V,n c 2b v,n “ ar V Sy, n 

b = — a±i 

v > v 4 8ln 9v 

b v , —«**&, 

,n (cos cpn— cos -cpv) [ 2(m+l) 


®V,v 


s v,n 


®V,n 


for 


_ _1 

2 (M+l) 


LA(v ; o)(fn l 9 ), + IySXy^+l). + Y L A (v,u) 

2 L—i 

t-i=l 


f n,U 


ar v = J3- = 

r cy chord at span station V 

L A (v,p), for a. straight wing, see equation (Ck) . 

L A (v,u), for a swept wing, see equation (C26) . 
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m 

1 1X1 •* 


n p-i 


^1=1 


cPn coa M-i 


■where again 


<Pn 


= ...nfl, 
n+1 


cpv 


vit 

m+1 


% 



The computations required, in the preceding group of equations 
may be simplified if such values as by v "bv,n f° r variolas m* e, 
and fn^iJ. for various m‘s and M*s are tabulated. Then a solution 
for any wing consists of a substitution of wing geometry into the 
L^C v ,n ) function, equation (C4) or equation (C26), and a combination 
of the tabulated coefficients to obtain m simultaneous equations 
with m unknowns Gq . 

The computations for a symmetrically loaded wing may be still, 
further reduced by an alteration to the preceding equations and 
coefficients. For a symmetrical wing with or without camber and 
twist, the distribution of local angle of attack is symmetrical 
about the center of the wing or 


then 


and 


a v + “m+i-v 

a, v 

Oy = CCj2-fl— V 

Gy = Gm+i— V 

v goes from 1 to 

n goes from 1 to m 
h goes from 1 to M 
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She summation terms can be written such that n from 1 to m 
becomes n + m+l-n from 1 to and. p from 1 to M becomes 

l-i + M+l-p. from 1 to Mi. (For the cases n = and p = ^S±i 

2 2 2 

the coefficients are expressed as their former values as will be 
seen directly.) 


Wow if equation (C20) is expressed as 


c. : 

= 0>*v,n + B*V,M+i-p) Gv- ' K,n + T>*v,m+i-n) G n 

n=l 


where 


n 4 

.2 

V = 1,2, B±l. 

2 


Then 


°v 


B\ Gy 


m±l 

2 

- Y, B * v >n 

n=l 


(C23) 


where 


B*v = B*v,p + v = 1,2 

B*v,n = b*v, n + b*^^ v « 1,2 


m+1 ,, 1 M+l 


ffi+1, n 1 Eli 
2 2 


(C29) 


These coefficients may also be expressed as 

B^V = 2by J v + Sv, V 
B*v,n = 2Bv,n - axv Sv ; n 


1 


(C30) 



•* 


^® re B v,n is limited in that B v#n » l» Vjn + fc %a+JUn for B v#n * \,n for n « ^ 

To find g v#n consider the egression for g VjT1 


gv = 1 

Sv ' n 2 (M+ 1 ) a 


M 

]T i a (v,m+i) f njlI 


(i=x 


The summation term of 8v, n can written for n 4 


M-l 

2 


^ L A^ v ^^n,ix = ^ ^ v ; n ) f n , ^ ) f n j ^ ^ I ■*■ ) f m+i— n j n ^ A^ J ^ + l~i i ) f m+i— ti , M+ i-jj. ^ 


|i=l 


<L-> 

H=1 


LA(v ^ )f m + x~n >M±L 


for n = S±1 
2 


*fci 


Ba( v ^m) fn,|i a / fiA( v ^u)f nj jji + IaCVjM+Ih 1 ) fn,M+i-ii^ + r X a( v f a M+i . 


fi-1 


H=1 


2 J 


(C31) 


VJl 

VJI 
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■where as before 


m 


f n,H “ 5+T 2, tAl Sln 111 m?X 008 8*1 

Hi=l 

examination of equation (C32) will indicate that for n = 


(C32) 


I *( v ’ I 2 L ) f n^ j 


= 0 


for n 4 

2 


have values only for even m.i, hut for even [i 1} f njM , = - fm+i-^ 1,11 
so that these terms equal zero. 

-y. . 

Then equation (C3l) becomes, for n —rr; 

M— 1 

V~; ^ 

^ l a(v,h) f njM . = 2_, r^ f n 3 u + f m+i-n,|i) La( v ^u) 


A., 

V-=l 


(JL=1 

+ (fn^M+i-ji + fm+i— n,M+x-ti) L A (v,M+l-fJ.) ] 


and for n = ®?p=> 
M 


M-l 

2 


^ L A (v,a) f nj n - [(^njl-i) Ia(v,u) + (^n^M+i-fa) £a(v,M+ 1 ~).i)] (C 33) 

m. 


Further examination of equation (C32) will show that for only the 
odd Hi 


f n,ii - — ^n^m+i-ix 



However, the second summation in equation (C33) has n = ^5^ in which instance the even nj. 
terms in equation (C32) vanish, thus the second summation oecomos 


ifci 

2 

l 

n=i 


f m+l a [Ijd(v,ji) — Ia(v,M+ 1 -p)] 


Tho first summation in equation (C 33 ) has terms of f n „ = - f a+lHa „ and fn u+l-n « - f m +i~n,M+i-u 

tr/an 1 » - am +V1 A r* AA'P^'l A*? AVI + A A^P T . / A I II ^ T . t A I I I I I ^ nil <PnM nirnM ^ ^ ‘ 


for even Pi 
In addition 


7 \ r “Jr * iut -‘-‘''t.i ~ - ’i • r'jrrj."* 

or the coofficients of %(V,|i) and L^(v,M+l-{i) vanish ror oven pi. 

f n,M+i-iJ. +:f 'm+i-n,M+i'ii for ^1 tho summation is 


fn,n +f m+i-n,iJ. = 


M-l 

~2 

( f n,p%+iHa,|J [ l a( v ^) ~ L A (v,M+l-ii)] 


where fn,|i + %+t-n,p, is obtained from (C 32 ) for the odd terms of Pi. 



lastly. 


(fn,o)l'A(v j0 )+(fn,M+i)I-A(v J ^ ) (fn,o+4n -i-n,o)l , A( V 5 c )'Kl , n,M+i +; ^KH-J.- j n J M+i)l , A^ v jM+1) . m+1 

2 “ 2 ° r n f 2 


B {-iuo +f m+iH i,o) M Vj °). La ^ V jM+1 ^ for n = ttti 
2 2 

where f n ,o +f m+i--n,o I s taken for only odd | ± x . And, for oven Pi f n , 0 - %+i-n,o. So that 
a factor f n#(i may he expressed as 

F n ,ix a f n ^ for n « and odd Pi 

f n,p ~ f n,P +f m+i--n,p for n ^ and odd Pi 
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and 


Examining equation (C32) 


fn,|i = for n4~ f H=0 and odd Hi 

fn,H = for n = H=0 and odd Hi 

tk± 

2_ 

^ VjT1 = 2(M+l) ) ^n,H ^A^jh) — La(v,M+1-h)] 

H=0 


f n.u ~ sin ^1 - J3 f 
m+1 


)& 


(C3h) 


fjn+i-n,H ** sin Hi 


(m+l~n)jt 

m+1 


sin 



sin Hi* cos - sin cos Hi* 
m+1 m+1 


or ’when Hi is odd. 


f m+i— n j H ru 


sin 


diDH 

m+1 


Thus 


fn,H ” ^n+i— n,H 
M— 1 

ev,n = V f n ,H Cla(v,u) - La(v,M+1-h)I 

H=d 


Thon, as before 
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where finally 

fn 3 p = for n = and- odd pi 

fn 3 p. = 2fn,[i for n ^ and odd pi 

fn,M- = fn,p. for n ^ ®j=pS y = 0, and odd m-i 

fn 3 p. = for n = p. = 0 3 and odd (Jt x 




) (035) 


and 

m 

fn.n = ) Pi sin Pi 22— cos pi ti2— 

n ’^ m+1 L, 1 1 m+1 M 

Pl=l 

A further simplification to g\, n can also he effected. From 
the binomial theorem: 


-a*iS=-4S*fe5S-*S + - 

! 1 + i i) _ h (i) + A 0D ~ A 


= a 


for a >b 


-] 

(036) 


and 


^ aS + b2 = b+ lh^“lt^ + ih 5 -'- 


’* [ a + I (tg) 


5 <5^) + A Vbfef 


jaffCs}^ + 


for a<h 


(C37) 


Now letting 


T(.t) = VZH! 

8 > 
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then for t = & 


T(t) = 1 + J t2 - 1 & + t? ”l2g t8 fort<1 


or 


T(t) . t +!('!') -iflY 3 + i- 

2 8 \t/ 1 6 


(t) = ~zk ft) 7 for 


Equation (C34) may be written 

M— 1 
2 

Sv,|i “ 2’(M+l) f n,|i 

H=0 

where 

£Ea(v,|i) = La( v ,u) - LA(v,M+l-y-) 
This can also be written 


= LaCt^) - LA(q,-n) 


so that from equation ( 025 ) 

AL A (v li) = ^l-^v ( n-^) tan_A ] g + [ary ( g-g ) 3 g 
* arv(q-Ti) 

1 J [l+ar v (i)— Tf) tan A] a + [ ar v (r)-t-T])] 2 

l+2ar v q tan A ar v (q+rf) 

g tan_A * — _ 

l+2arvq tan A */ [1+arvq tan A] + [ajcvn ] 

25 

ar v (i]2- q2) 


t>l 

( 038 ) 

(039) 

(C40) 

(c4i) 


(C42) 
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Then 


AL A (v,h) - T( t x )- 


l+2aryr|tan A 
- T(t 3 > tan A - 


T(tg) +T(ta) tan aJ 

, / ; 


2 a. 


ary 


) 


(C43) 


where 


ti 

tg 

ts 


atv(n-n) 


zsr- + tan Al 


1 — „ -t- Ifc& tan A 

ary (q+T|) t.+ti j 

ary T) 


l+ar v tan A 


and T(ti) takes the sign of ( tj— T f) 

T(t 3 ) takes the sign of t 3 

Several of the functions of equation (cit-3) are independent of wing 
parameters and may "be tabulated for various tj and tf. 


Let |3v 


1 

az\)* 





K, = K = ^ 

* r]£— ri£ 5 ri+q 


•where, as before. 
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Then 


AL A (v,n) = T(t x ) - 


1+2 tan A 


[ T(t 2 ) + T(t 3 ) tan A] 


Pv 2: S 

T(tg) tan A - P v E 4 


where 


and 


av 

B* v 

B* 


V,n 
Bv,n 

*>v,n 


ti — jp v Ei + tan Aj 

t 2 = |Pv E 2 ■+■ Eg tan Aj 

t 3 = — 

Pv Kg + tan A 

T(ti) takes the sign of Ki 

T(t 3 ) takes the sign of t a 

In summary, for the symmetrically loaded wing, 
m+1 
2 

= B*v Gv ~ y B* v ,n Gn, V = 1,2 . . . 
m=l 

= 213-Vj v + Gv, v 
= ^Vjn “ ar v s v,n 


r h v, n 

+ h^ 

;,m+i— n f° r 

n 4 m ^ L 


for 

n = B±i 
n 2 


jp£ L 




4 sin 

*9v 



r 

sin 

CPn 


_ (cos 

<Pn ' 

- cos cp v ) 2 

2m+l J 


(c44) 
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g v n is taken frcsn equation (C39); ^LAC'Vjt 1 ) is taken from 
(C&i-) and f njli j_ s frcsn equation (C35) * 

Limitation of tiie Series 

The number of coefficients G n required for accuracy depends 
upon how rapidly the series equation £ 3 ) converges. Weissinger 
used m equal to 7j 15 j and 31 in his investigation and concluded 
that the results obtained with m equal to 7 were nearly as 
accurate as those with m equal to 15 or 31- For this reason m 
equal to 7 > or four coefficients have been used in all of the applica- 
tions of this method presented herein Hie number of terms required 
in the interpolation function f n ^ must also be established. Again 
Weissinger used M equal to 7, 1$, and 31 and found that results 
with M equal to 7 proved as satisfactory as those with M equal 
to 15 or 31 - Lastly, it should be noted that equation (3) cannot 
satisfactorily approximate a. curve containing discontinuities; 
however, a modification which will enable it to do so has been 
developed by Multhopp (reference 10) . 


Solution for Additional Loading 

Since in a solution for additional loading the wing is 
considered a flat plate and all angles of attack cty are equal to a, 
equation (C28) may he modified to 

m+1 

1 - B *v ~ - £ B *v,n I 1 Cc45) 

n=l 

Evaluation of this equation at the several stations v produces a 
set of equations containing the unknown circulations Gn which can 
•then be solved simultaneously to obtain the values of these 
circulations . 

Substitution of the values so obtained into the following 
expressions, results in the values indicated: 



n=l 


(C46) 
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or for a symmetric wing, 

m-1 

4CL = SAE^^IZ +2 f‘ 2a 

da, m+1 \ a, / , ~~ 

n=l 


f 1 si ^ % 


(C 47 ) 


ci ? G v 1 

— = 2 ary — 


a dC^/dcc, 


- I - C - v ' = 2AB ~ — 

^L°&v a, dCL/da 


(C48) 

( Ck9 ) 


and from reference 11 for m = 7 


^cp 


0,3$24Gi + Q.5Q3G g + 0.3^10-3 + Q.0405&4 
G.3827G1 + 0.7071G-2 + 0.9239G3 + O.5OCG.1 


(C50) 


Solution for Basic Loading 

The basic loading on a wing with camber and/or twist can be 
determined in a manner exactly parallel to those of Falkner and 
Mutterperl. An arbitrary angle c% is selected for the root section 
and the values a,iocal> measured from it. If these values are then 
substituted in equation (C28), and if the resulting equations are 
solved simultaneously, values of G n will be obtained, which when 
Inserted into the following expression will give the correct lift 
coefficient for the wing at this attitude. 

Efci 

C L = ^ +2 

" m+1 ^ 2 

n=l 

If these values of G n are also substituted into the following 
equation, an expression for the total loading curve will result. 


1 


Gn sin cpn 


(C51) 


ci cy = 2b Gy 


(C52) 
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Of the values of Gq obtained frcaa calculations for the addi- 
tional loading aro thon substituted into equation (C 52 ) , ordinates 
of the additional loading curvo will be obtained,, which when 
subtracted from those of the total loading curve will givo those 
of the desired basic loading curvo . 


Correction of ci a 


As in the other methods, the error introduced by the assumption 
that the section lift-curve slope is in all instances 2it can be 
readily corroctod. Tho correction is accomplished by modifying the 


specific circulation ordinates Gn by the ratio 

specific values of this function at span station 
determined if this function varies along tho span. 


actual ci a 

iit 

n must be 


whore 


Computing Instructions 

The AL^(V,ii) functions are determined on form c(l) for a 
swept wing and form C(2) for a straight wing (Ac. 25 c = °) • In 
both cases the coefficients Ki, K2, K3, Sit, and K 5 aro obtained 
from tablo CX. in tho values of T(ti), T(t 2 ), T(t 3 ) aro obtained 
by entering chart Cl with the values of t x , t 2 , t 3 from columns 
(10), (ll), and (12} form C(l) . Form C ( 3 ) contains thoso computa- 
tions which result in the osprossions containing tho 

unknowns G^. The values gv,n column 9 of form C ( 3 ) one obtained 
as follows: Consider tho valuos in column J+ of this form as four 

groups of four members with tho groups identified as 

A when n = 1 

B when n = 2 

C when n = 3 

D when n = h 

Similarly, tho valuos in column 5 of this form can bo considorod 
grouped as 
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1 when v « 1. 

2 when v = 2 

3 when v = 3 

4 when v = 4 

then 

Sx,X = 2 A X 1 
Si, a = 2 B X 1 
s = ^ C X 1 
& 2.^4 = S D X 1 
3a, 1 = 2 A X 2 
e SfZ - 2 B X 2, etc . 

The computing form for the simultaneous solution (reference 9) 
of the equation (C45) is given in forat C(4) . The equations are set 
up as follows: The first equation consists of the first four numbers 

In column 15 form C(3), the second equation the second four, and so 
forth; the first number in each group being the coefficient of Gj./a,- 
the second being the coefficient of G^/a and so forth. Simultaneous 
solution gives the G/a’s with tha corresponding span station. 
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Form C : ~ Computing form for Weissingers method 

( Underscored numbers are sample calc ufc i lOifS.) 
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TABLE Cl.- CONSTANT FACTOE ESQUIRED FOB THE 
CALCULATION OF £L A (v,[i) FCS m = M = 7 



M- 

Ki 

z 2 

Ks 

K 4 

K s 


S-AHT 

i/n+n 

l/r] 


tHT/ti+tT 

1 

0 

- 13.1406 

0.5198 

1.0824 

-13.6610 

-0.0396 


1 

00 

.5412 

1.0824 

CO 

0 


2 

4.6125 

.6131 

1.0824 

3.9992 

.1329 


3 

1.8477 

.7653 

1.0824 

1.0823 

.4142 

2 

0 

- 3.4141 

.5858 

1.4142 

- 4.0000 

-.1716 


1 

-4.6125 

.6131 

1.4142 

-5.2254 

-.1329 


2 

CO 

.7071 

1 . 4 i 42 

00 

0 


3 

3.0826 

.9176 

1.4142 

2.1650 

.2977 

3 

0 

-1.6200 

.7232 

2.6130 

-2.3432 

-.4464 


1 

- 1.8477 

.7653 

2.6130 

-2.6128 

— , 4 i 42 


2 

-3.0826 

.9176 

2.6130 

- 4.0002 

-.2977 


3 

“ 

1.3065 

2.6130 

00 

0 

4 

0 

-1.0000 

1.0000 

CO 

- 2.0000 

-1.0000 


1 

- 1.0824 

1.0824 

00 

- 2.1649 

-I.0000 


2 

- 1 . 4 i 42 

1 . 4 l 42 

ra 

- 2.8284 

-1.0000 


3 

-2.6130 

2.6130 

00 

-5.2260 

-1.0000 
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TABLE CXI.- INTERPOLATION FUNCTION 
AS CALCULATED WHE a * M « 7 


n 

V- 


1 

0 

2.613 

. 

1 

-1.414 


2 

-1.531 

■ 

3 

.414 
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0 

— 1.4l4 


1 

2.696 


2 

- 4.000 


3 

-1.531 

3 

0 

1.082 


l 

-2.414 


2 

3.696 


3 

— 1.4l4 
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- .500 


1 

1.082 


2 

-1.414 


3 

2.613 


TABLE CIII.- VALUES OF B v , n 
AS CALCULATED WITS m = 7 


V 

n 

2 Bv# 

l 

1 

0 


2 

3.8284 


3 

0 


4 

.2928 

2 

1 

2.0720 
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0 


3 

2.3888 


4 

0 

3 
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0 


2 

1,8284 


3 

0 


4 

1.7022 

4 

1 

.2242 


2 

0 


3 

3,1548 


* 

0 


TABLE CIV.- VALUES OF L VjV 
AS CALCULATED WHS m=? 


r ~v ~ 1 

nr 

2 b v . v 

1 

1 

2 

3 

4 

10.4524 

2 

1 

2 

4 

5.6568 

z z z z 

3 

1 

2 

3 

4 

4.3296 

4 

1 

2 

3 

4 

4.0000 
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